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Abstract 

We measure the horizontal (| cos(# 2 )| < 0.14 corresponding to a slant depth cut 14 kmwe) 
neutrino-induced muon flux (E^ > 1.8 GeV) in Soudan 2 to be 4.01 ± 0.50 ± 0.30 x 
10~ 13 cm _2 sr _1 s _1 . From the absence of horizontal muons with large energy loss, we set 
a limit on the flux of muon neutrinos from Active Galactic Nuclei. 



1 Introduction 

The Earth provides an effective shield against the muon component of cosmic ray 
showers for overburdens exceeding 14,000 meters-water-equivalent (mwe). Conse- 
quently, near the horizontal direction (| cos(6> z ) | < 0.14) at the site of the Soudan 2 
detector (vertical depth 2090 mwe), the muon flux is comprised almost entirely 
of initiated by charged current interactions of atmospheric muon neutrinos in 
the rock surrounding the detector. In addition, astrophysical neutrino sources may 
contribute to the (horizontal) muon flux. In particular, a high energy (> 5 TeV) con- 
tribution may arise as the result of neutrinos which originate within astrophysical 
sources such as Active Galactic Nuclei (AGNs)[ 1,2,3,4,5, 6,7]. 

Active galaxies are those with extreme variations in brightness corresponding to 
abnormal emission of large amounts of energy at optical and/or radio wavelengths. 



Preprint submitted to Elsevier Science 



7 February 2008 



Colliding and exploding galaxies also fall into this category. The nuclei of these 
galaxies may produce both neutral and charged particles. Most of the particles 
either decay or are absorbed before they can escape the galaxy, leaving only photons 
and neutrinos to be detected on Earth. Several models predict substantial high- 
energy neutrino production from AGNs[l,2, 3,4,5, 6,7], however there is a theoret- 
ical limit which conflicts with many of these models [8]. AGNs may collectively 
give rise to a diffuse neutrino flux, or they may be detectable as individual point 
sources. Gamma Ray Burst sources have also been proposed as another possible 
diffuse source of high energy neutrinos[9]. 

Charged current interactions of multi-TeV muon neutrinos (y^ + v^) in the rock 
surrounding the Soudan 2 detector can be expected to give rise to TeV muons at the 
detector. Muons with these large energies will very frequently exhibit catastrophic 
energy loss by bremsstrahlung or pair production as they traverse the tracking 
calorimeter. This is highly improbable for muons (< >~ 20 GeV at the 
detector) initiated by atmospheric neutrino reactions in the surrounding rock. Thus 
it is possible to measure the rate of high energy AGN neutrinos. 

In this paper we report a measurement of the horizontal muon flux of 4.01 ± 0.50 ± 
0.30 x 10~ 13 cm _2 sr _1 s _1 and a negative search for high energy muons from AGNs. 



2 The Soudan 2 Detector 

The detector is located 
approximately 700 m 
below the Earth's sur- 
face in an historic iron 
mine in Soudan, Min- 
nesota, USA. A pho- 
tograph is shown is 
Figure 1. 

The central detector 
is an iron sampling 
calorimeter using drift 
tubes filled with 85% 
Ar - 15% C0 2 gas 
as the active medium. 
The detector is made of 224 modules with dimensions of 2.7 x 1 x 1.1 m 3 . Each 
module has a mass of 4.3 t. The modules are stacked 2 high for a total height of 
5.4 m, and 8 wide for a width of 8 m. The average density is 1.6 g/cm 3 . 

The drift tubes are oriented along the North-South direction and are operated in 
proportional mode. Each end of the tubes is read out by an anode wire and a cathode 
pad which establish the x (East- West) and y coordinates of the hits. Off line, hits 




Fig. 1. The Soudan 2 iron tracking calorimeter is a 8 m x 
5 m x 16 m assembly of standard-unit 4.3 t modules. Each 
module contains corrugated Fe sheets interleaved with drift 
tubes and stacked into a fine-grained hexagonal lattice. 
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are matched in each half tube, using pulse height to match the cathode and anode 
pulses. The z (North-South) position of the hit is determined by correlating the drift 
times. 



The ceiling, floor, and walls of the detector hall are covered with an active shield 
consisting of a two layer aluminum proportional tube array. The shield array re- 
solves charged particle crossings ("hits") to within four microseconds and within 
rectangular surfaces of 10 cm x 5 m. For most Soudan 2 analyses, the shield is used 
to veto events in the main detector which were initiated by particles interacting in 
the surrounding rock. Use of coincident in-time shield hits in the present analysis 
is described in Section 4. 

The trigger for recording events requires a multiplicity of seven anode pulses or 
eight cathode pulses in blocks of 16 contiguous readout channels. For long muon 
tracks this trigger is fully efficient. 

A near horizontal muon data event which is (presumably) neutrino-induced is shown 
in Figure 2. This example is the longest horizontal muon track that was found, with 
a length of 1464 cm and 149 tube crossings used in the fit. The zenith (azimuth) 
angle of this track is 84° (23°). Two shield hits can be seen in the top and front 
views. 

The construction and performance of the tracking calorimeter are extensively de- 
scribed in Reference [10]. Of particular interest to this analysis, the single track an- 
gular resolution has been measured using the moon shadow[l 1] to be (0.3° x 0.3°) 
and the electromagnetic energy resolution determined by counting tube crossings 
has been measured in the ISIS test beam[12] to be 

where E is given in GeV and the second term reflects the saturation of the number 
of hits. 



3 Horizontal Muons at Soudan 2 

Figure 3 shows the vertical muon intensity underground as a function of depth in 
standard rock, as compiled by Crouch[13] and subsequently updated by the Particle 
Data Group[14]. Two distinct components are apparent in the fitted curve, which 
consists of a double exponential plus a constant term. The atmospheric muon rate 
is observed to fall sharply with slant depth, so that the flatter spectrum of neutrino- 
induced muons only becomes visible for slant depths greater than 14 kmwe. For 
angles 9 Z with respect to the vertical (zenith angle), the slant depth increases ap- 
proximately as sec(0 2 ), but depends in detail on both the surface terrain and the 
rock density and composition. 
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TOP VIEW 



Soudon2 Data 
Run 67014 Event 761 
1 0-Feb-1996 07:06:31.87 
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FRONT VIEW 



Fig. 2. Three views of a neutrino-induced horizontal muon traversing the Soudan 2 detector. 
The muon traverses eight calorimeter modules. The proportional tube shield registers its 
entrance and exit points on the cavern walls (see TOP and FRONT views). Dimensions are 
in cm, and individual module outlines are shown. All three views use the same scale so that 
a correct aspect ratio is maintained. 

Based on the vertical muon spectrum and the calculation of the average slant depth 
at Soudan 2 as shown in Figure 3, we can select a sample of neutrino induced 
muons by making a cut on the zenith angle near 82°. Since the event timing res- 
olution is about 100 ns, it is not possible to ascertain whether a muon is upward- 
or downward-going. There is thus a two-fold ambiguity in the direction for each 
measured track and 9 Z is defined to be < 90°. 

Calculation of the slant depth requires knowledge of both the rock density and the 
surface terrain. Differences between standard rock (p ~ 2.62 g/cm 3 ) and Soudan 
rock (3.00 g/cm 3 > p > 2.74 g/cm 3 ) must be taken into account. The topology 
above and surrounding the Soudan site is shown in Figure 4[15]. The geology of 
the region is such that the iron deposits are almost vertical, so that the most vertical 
zenith angles correspond to the highest average density. The majority of the rock is 
greenstone composed of silicon dioxide (Si0 2 ) and aluminum oxide (Al 2 2 ).[16] 
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DEPTH (m.w.e.) 

Fig. 3. The vertical muon intensity as a function of slant depth, from the compilation of M.F. 
Crouch[13,14]. Arrows indicate the average depth at different zenith angles (in degrees) at 
the Soudan 2 detector. 

Based on Soudan rock effective parameters, a depth of 14 kmwe in standard rock 
corresponds to 13.2 kmwe in Soudan rock.[17] 

Due to the uneven surface terrain, the slant depth was calculated for each track 
individually. Observed fluxes of cosmic ray induced muons with 6 Z < 60° have 
been used to determine the zenith angle dependence of the rock density[18]. By 
considering the azimuthal variation in the slant depth, we infer that systematic 
uncertainties in the atmospheric muon flux due to incomplete knowledge of the 
rock density are on the order of 10%. 

Upon inclusion of the terrain information, a more refined estimation of the slant 
depth as a function of zenith angle is obtained, as indicated in Figure 5. The solid 
curve in that Figure was made assuming a spherical Earth. The dots which appear 
each one degree in zenith (two degrees in the inset) and are averaged over ten 
degrees in azimuth reflect the effect of topology variations in the overburden. 



4 Reconstruction of Neutrino-Induced Muon Tracks 

There were approximately 10 8 triggered events in the Soudan 2 detector in the data 
set considered for this analysis. About half of these were due to cosmic ray muons, 
and half due to radioactivity in different parts of the detector which satisfied the 
trigger due to multiplexing. These events were processed using a pattern recogni- 
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1 Degree Digital Elevation Model (1 :250,000 Scale) near Soudan MN 




-93.0 -92.8 -92.6 -92.4 -92.2 -92.0 -91.8 -91.6 -91.4 -91.2 -91.0 

01/04/97 16 44 Longitude (degrees) 

Fig. 4. USGS 1 degree digital elevation model extending to Lake Superior. The zenith 
angles 78 < 9 Z < 90 are projected onto the surface. The distance between Soudan and 
Two Harbors, MN, is 97 km. 





Data Selection Criteria 


1 


Track Status (extra-fiducial detector hits and in-time shield hits) 


2 


Track Length Cut (> 175 cm) 


3 


Azimuthal Cut (> 8° from quadrant axis) 


4 


Zenith Cut (> 78°) 


5 


Anode Width 



Table 1 

Selections used to define the horizontal muon sample. 

tion filter code which identifies muon tracks and discriminates against unphysical 
background from noise[17]. Event processing took place in three stages. Runs were 
processed by the mine computer soon after the runs were complete, and sorted into 
a number of output streams, one of which was a file consisting of events which 
included a track with a zenith angle greater than 60 degrees. Those files were then 
subjected to a filter program described below which identified horizontal muon 
candidates with a zenith angle greater than 78 degrees. There were 4458 events 
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from the output of that filter which were scanned. 

The filter program used five cuts which were applied to define a "horizontal muon" 
sample as listed in Table 1 . Tracks were required to have end points consistent with 
a particle entering and leaving the detector. Acceptable track ends were required 
to 1) lie outside the fiducial volume which is 50 cm from the edge of the main 
detector, or 2) lie on a module boundary inside of the main detector and 3) point to 
a coincident shield hit or 4) point to a portion of the shield where a shield module 
did not exist. The azimuth cut rejected noisy events reconstructed parallel to the 
axes of the detector. The zenith cut defined the sample of horizontal muons. The 
track length cut minimized short tracks such as pions that originated from nuclear 
interactions within the rock nearby and also low energy cosmic ray muons which 
had undergone large multiple Coulomb scattering. 

The anode- width cut re- 
jected events that were 
primarily along one an- 
ode wire. These were 
actually along the north- 
south direction, but were 
often mis-reconstructed. 
The efficiency of these 
cuts for Monte Carlo 
muons generated isotrop 
ically with zenith angles 
exceeding 82 degrees is 
0.56 ± 0.01 as given in 
Table 2. Note that both 
the trigger and recon- 
struction efficiency are 
lower near zenith angles 
of 90 degrees. 

The pattern recognition 
software found single 
horizontal muon tracks, but it also kept a variety of background events. Most 
background events contained short horizontal tracks accompanied by other activity 
in the detector. Physicist scanning of all tracks greater than 78° was done to 
verify the track fits and to eliminate obvious backgrounds. Common backgrounds 
were due to certain patterns of electronic noise, vertical muons undergoing large 
radiative stochastic losses, multiple muons, and neutrino interactions within the 
detector. Tracks with visible multiple scattering (greater than 2 degrees) or with 
in-time shield hits that were not due to a horizontal muon were also rejected. All 
events were subjected to two independent scans and residual discrepancies were 
rescanned. The uncorrelated scan inefficiency as determined from the double scan 




85 90 
Zenith (degrees) 



Fig. 5. Estimation of slant depth versus zenith angle for a 
spherical Earth (curve) and upon inclusion of terrain and 
density information for the Soudan site (dots). The slant 
depth is shown in units of meters (left) and overburden 
(right). 
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was less than 1%. The Monte Carlo which was used for the efficiency calculation 
only generated high energy (100 GeV) muons. This fails to take into account 
neutrino induced muons with lower energies in the detector which might fail the 
multiple scattering cut applied during scanning. An independent calculation that 
took into account the neutrino induced muon energy distribution at the detector 
gave the result that 5% ± 5% of the events should undergo multiple scattering by 
more than 2 degrees in the detector, the criterion used during scanning. We therefore 
used a scan efficiency of 95% for the atmospheric neutrino flux calculation. Note 
that this inefficiency does not apply to the AGN v search in the next section. 



Cut 




Generate 


2000 


Matched Track 


1920 


Track Status 


1818 


Track Length 


1685 


Azimuth 


1411 


Zenith 


1378 


Anode-width 


1111 


£ 


55.6% 



Table 2 

Survival rates of 2000 Monte Carlo muons incident nearly horizontal on the Soudan 2 
detector. 

The track-length cut of 1.75 m (2.2 hadronic interaction lengths) minimized the 
number of short tracks which arose from processes which were not neutrino-related. 
The latter backgrounds included pion tracks originating from deep inelastic muon 
scattering within nearby rock, or else high energy cosmic ray muons from vertical 
directions which experienced large-angle multiple scattering. This length corre- 
sponds to a muon energy threshold of 400 MeV, but the restriction on multiple 
scattering raised this to an effective threshold of 1.8 GeV. 

The number of accepted events after scanning, with 9 Z > 78°, was 1237 events. 
The slant depth cut of 14 kmwe reduced this sample to 65 events. The track length, 
azimuth and zenith angle distributions of all 1237 events are shown in Figure 6. 
Figure 7 shows the (9 Z , 0) distribution of these events. The depletions at regular 
intervals on azimuthal projections of these plots are' due to the azimuth cuts. The 
wavy-line contour on Figure 7 shows the 14 kmwe slant depth cut, which is de- 
signed to separate the sample into atmospheric muon and neutrino-induced muon 
candidates. The event with the largest zenith angle is shown in Figure 8. 

To estimate the amount of background from atmospheric muons we integrate under 
the intensity-vs. -slant-depth curve of Figure 3, but without the constant (neutrino- 
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Fig. 6. The reconstructed zenith, azimuth and track-length distributions for through-going 
muons which were selected by the filter software to have a zenith in excess of 78° and 
passed the scan criteria. The shape of the track-length distribution reflects the shape of the 
detector for an isotropic muon flux. 



induced muon) term, for slant depths greater than 14 kmwe. This procedure yields 
a background estimate of 0.1 l^i events in our horizontal muon sample, where the 
error reflects the uncertainty in the rock density. 

Muons which multiple scatter in the Earth above the Soudan 2 detector are a poten- 
tial source of additional background. Since the overburden above the Soudan 2 
detector is so great, it is not appropriate to apply the usual equation for small 
angle multiple scattering. Instead, a Monte Carlo calculation was performed which 
propagated muons through the Earth including multiple scattering and calculated 
the energy loss in steps. With a sample starting with 60 million muons that had 
9 Z > 60 degrees at the surface of the Earth, the same number that were recorded 
during our exposure, no additional muons were generated with 9 Z > 82 degrees 
at the detector. The finite statistics of the Monte Carlo introduced an additional 
uncertainty which corresponds to < 0.5 event background. We conclude that the 
expected background is less than one event. 

The neutrino-induced muon flux, <fr^, can be expressed as 
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Fig. 7. Angular distribution of horizontal muon candidates. The contour line near 9 Z = 82 
represents the 14 kmwe contour. 



N 

where A^ t is the number of observed neutrino induced muons, t is the time exposure 
of the detector in seconds, A is the effective area of the detector in cm 2 , f2 is the 
solid angle subtended by the detector in steradians, and e is the detection efficiency 
for the muons in the detector. Both A and are calculated assuming a uniform 
acceptance to the right of the contour in Figure 7. 



The analyzed data spanned the period from April 14, 1992 to April 24, 2002. The 
exposure time was calculated from the start and end times of every processed data 
run. The calculated exposure for this analysis, including corrections for the detector 
duty cycle and electronics dead time during data taking, is 2.00 x 10 8 s. To estimate 
the effective area for this analysis, Monte Carlo muons were uniformly generated 
and the effective area for each track calculated. For 9 Z > 82° the average effective 
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TOP VI EN 



Soudan2 Data 
Run 58027 Event 337 
23-Dec-1994 00:36:17.03 



Fig. 8. Our largest zenith angle event. 9 Z = 89.9° and the length is 731 cm. The scales are 
labeled in cm. 

area was 86.7 ± 0.3 m 2 . The solid angle was calculated from the acceptance region 
in Figure 7 to be 1.77 sr. The trigger and reconstruction efficiency e = 0.53 (0.556 
from Table 2 x 0.95 from scanning.) 

Systematic errors arise with the uncertainty in the background of 0.7 events (1%), 
uncertainties in the scan efficiency (5.6%) and uncertainty in the energy distribution 
as it affects the efficiency of the 2° multiple scattering cut (5%). Adding these 
errors in quadrature, we assign an overall systematic uncertainty to the neutrino 
flux calculation of 7.6%. 

The resulting neutrino-induced muon flux for a muon energy threshold of 1.8 GeV 
is 

= 4.01 ± 0.50 (stat) ± 0.30 (sys) x 10~ 13 cm^sv^s- 1 (3) 

Upward-going atmospheric neutrino-induced muons have been previously mea- 
sured by Baksan, 1MB, MACRO, Kamiokande, and Super-Kamiokande[ 19,20,2 1,22,23]. 
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However, these experiments have not reported neutrino-induced muon fluxes above 
the horizon from those regions of solid angle where the overburden might be suf- 
ficient to separate atmospheric muons, perhaps due to uncertainties in the surface 
topology and rock density. The SNO experiment is deep enough that it should be 
able to see the neutrino-induced muon flux above the horizon[24]. The flux near 
the horizon has also been measured by LVD[25] and Frejus[26,27]. The present 
measurement has aspects in common with that of Frejus, albeit with increased 
statistics. LVD, which has a much less uniform overburden, did not separate their 
neutrino-induced signal from background. The measured fluxes closest to 9 Z = 90° 
are tabulated in Table 3 along with our estimate of the muon energy threshold used 
for each analysis. 



Experiment 


Flux 


E ■ 

i-'min 


cos 9 Z 


Depth 




(KT^cm^sr-V 1 ) 


(GeV) 




MWE 


Baksan[19] 


4.04 ±2.01 


1 


-0.1 to 0.0 


850 


Frejus[26] 


3.67 ±0.66 


0.3 


-0.18 to 0.18 


4710 


LVD [25] 


8.3 ±2.6 


1.0 


-0.1 to 0.1 


3000 


1MB [20] 


5.66 ± 0.95 


1.8 


-0.14 to 0.0 


1570 


Kamiokande[22] 


2.84 ±0.53 


1.7 


-0.1 to 0.0 


2700 


Super-K[23] 


3.45 ± 0.33 


1.6 


-0.1 to 0.0 


2700 


MACRO [21] 


7.4 ±2.8 


0.4 


-0.1 to 0.0 


3150 


Soudan 2 


4.01 ±0.58 


1.8 


-0.14 to 0.14 


2090 



Table 3 

Comparison of near-horizontal neutrino-induced muon fluxes measured by several 
experiments. Zenith angle ranges and energy thresholds are approximate. Several 
experiments report fluxes for restricted portions of azimuth angle to reduce backgrounds. 
The reported depths are minimum overburdens, except for Frejus and Baksan which report 
effective depths. 



Super-Kamiokande[28], and Kamiokande[29], MACRO[30], and Soudan 2[31] as 
well have analyzed atmospheric neutrino events and concluded that neutrino oscil- 
lations modify the zenith angle distribution of the flavor component of the atmo- 
spheric neutrino flux. For a neutrino energy of 40 GeV and assuming Am 2 ~ 3.5 x 
10~ 3 eV 2 and maximal (#23) mixing, the probability of oscillation for a from 82° 
(< L > ~ 130 km) is 2. 1 x 10~ 4 while for a from 98° (< L > ~ 2000 km) it 
is 0.04. Our modest statistics and 90° — 9 Z ambiguity make an oscillation analysis 
of these events impractical. 
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5 AGN Neutrino Search 



Several models for neutrino production in Active Galactic Nuclei predict large and 
potentially measurable fluxes of very high energy neutrinos[l,2]. These neutrinos 
would then produce high energy muons with energies from a few TeV up to 100s of 
TeV. In this energy range the dominant energy loss process for the muons in iron is 
electron pair production, followed by bremsstrahlung[32]. Both processes produce 
large electromagnetic showers in the Soudan 2 detector which are easily detected 
and measured. 



To understand the response of the tracking calorimeter to these TeV muons, a sim- 
ulation was performed. High energy muons were propagated through the Soudan 2 
detector and the amount of energy they deposited was calculated in a GEANT-based 
Monte Carlo[33]. Three different muon energies were studied: 5, 20 and 100 TeV. 
Figure 9 shows the fraction of muons which lost at least a specific amount of energy 
versus that energy loss for three different muon energies. The results obtained from 
the simulation agree with an analytic calculation carried out independently [34]. 
While 60% of the 5 TeV muons lose 5 GeV or more, 91% of the 20 TeV muons 
and 99% of the 100 TeV muons lose at least 5 GeV in the detector. 



The muons identified for the horizontal muon flux measurement were subjected to 
a predetermined cut of 5 GeV on the amount of energy loss they experienced in the 
detector. None of the 65 events had visible radiated energy loss greater than this cut. 
The largest zenith angle event with 
substantial visible energy loss is g 1 

o 

shown in Figure 10. This event i 09 
has a horizontal muon with a | o.s 
zenith angle of 80.6° and an | 07 
azimuthal angle of 133° which 06 
entered the top of the detector 

L 0.5 

and exits the north wall. Parallel 
to the muon at the top is an 
entering electromagnetic shower, 
presumably from a bremsstrahlung 02 
or pair production process along 01 
the muon path in the rock above the o 
detector. From Figure 7, the zenith 
angle (80.6°) tags the event as a 
high energy cosmic ray muon and 
not a horizontal muon candidate. 




120 140 160 

Energy (GeV) 



Fig. 9. Calculation of non-ionization energy 
loss for muons in Soudan 2 
Visible energy loss from a cosmic ray muon is not unexpected, because for 



atmospheric muons at the detector near zenith angle 80°, < > 



380 GeV 



while for atmospheric neutrino induced muons, < E^> ~ 20 GeV. The energy 
of the shower is reconstructed to be 2.2 GeV. There are no such events among the 
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65 horizontal muon candidates. 



Models for neutrinos from AGNs predict a v flux as a function of energy. For each 
assumed functional form of the predicted neutrino energy spectrum, an observation 
or limit on high energy muons can be used to normalize or limit the neutrino flux. 
It is less meaningful to plot a limit on a neutrino energy distribution for an arbitrary 
functional form. Therefore we choose to show our limit on the neutrino-induced 
muon energy distribution. For each muon energy, our non-observation of events 
with energy loss corresponds to an upper limit on the muon flux. 

To provide a comparison of our search with other searches reported in the literature, 
we show in Table 4 the exposures (AxtxVlxe). Our exposure is larger than that of 
the Frejus experiment, but lower than MACRO or AMANDA. Since the latter two 
experiments focus on upward-going muons, and the Frejus and Soudan 2 limits use 
horizontal muons, the limits are complementary for the exotic case involving the 
highest energy neutrinos for which the Earth's neutrino opacity is large (E u > 50 
TeV)[9]. 



Experiment 


Exposure (10 14 cm 2 s sr) 


Soudan 2 


1.71 


Frejus[27] 


1.60 


MACRO [35] 


3.40 


AMANDA[36] 


13.9 



Table 4 

Comparison of exposures (Product of running time, effective area, solid angle and 
efficiency) for AGN neutrino search experiments 

Based upon the observation of zero events, together with our efficiencies, we cal- 
culate 90% CL upper limits for the integral muon flux above three energies, which 
are shown in Table 5 and Figure 1 1 . Our limit is close to the highest predicted flux; 
however that flux is not ruled out. The limits are lower for higher energy cut-offs 
because of the higher probability of detecting a large energy loss. 



Soudan 2 Flux Limits for AGN Neutrinos 



Energy (TeV) 


High Energy (i Efficiency 


90% CL Limit (cm^sr^sec -1 ) 


5 


60% 


2.2 x 10~ 14 


20 


91% 


1.5 x 10" 14 


100 


99% 


1.4 x 10~ 14 



Table 5 

Efficiency for high energy muons to experience at least 5 GeV energy loss in the Soudan 2 
detector and the resulting AGN neutrino flux limits. 

Other experiments have also used their horizontal muon events to set limits on the 
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TOP VIEk 



Soudan2 Data 
Run 102038 Event 1030 
24-Oct- 1 999 09:44:44.68 



FRONT VIEW 




Fig. 10. Views of the muon event of largest zenith angle (80.6°) having a substantial 
radiative energy loss. The muon's zenith angle is too small to qualify as a neutrino-induced 
candidate. The vertical scale is magnified and all scales are labeled in cm. 

flux of the diffuse AGN neutrinos. The Frejus experiment used two methods to set 
a limit on the flux of high energy neutrinos [27]. One method extrapolated a fit to 
the atmospheric neutrino energy spectrum obtained previously [26]. Since that fit 
assumed that one spectral index fit all of their data, using it to deduce a limit for 
muons arising from a different energy distribution may bias that limit toward small 
values. The limit was quoted as d(j)/dE v (2.6 TeV) < 7.0 lO^GeV^cm-VW- 1 . 
With their second method, Frejus used the absence of muons with large energy loss 
to constrain the normalization of a number of specific models of the high energy 
neutrino flux[27]. 

The MACRO experiment[35] conducted an AGN search and found 2 candidate 
events with background of 1 . 1 expected. They set an upper limit on the muon flux 
from diffuse neutrinos of (1.7 ± 0.2) x 10 _14 cm _2 s _1 sr _1 . Their efficiency was 
highest for upward-going muons and lowest for horizontal muons. The AMANDA 
neutrino telescope, located in the ice near the South Pole, has set a limit based on 
an assumed diffuse E~ 2 spectrum of dN/dE v < lO^E^cm-V^r^GeV^pe] 
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As another test for an extraterrestrial component in the horizontal muon events, we 
have examined them for evidence of point sources. Figure 12 shows the directions 
of the 65 horizontal muons in galactic coordinates using an Aitoff projection of 
the galaxy. For each muon, the two possible directions are plotted. We observe that 
no clusters involving two or more muons within the detector's angular resolution 
(0.3° x 0.3°)[11] appear within the plot. 
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Fig. 11. The data points are our totally correlated limits from Table 5. The Curves a-f are 
from horizontal muon flux predictions from models by Szabo and Protheroe [2] and curve 
g from Stecker [5]. The shaded area represents the atmospheric neutrino flux. 



+90 
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Fig. 12. Aitoff Projection of the 65 horizontal muon candidates. Each event has two 
projections due to the ambiguity in direction. 
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6 Summary 

In this study we have isolated a sample of 65 neutrino-induced muon candidates, 
with an estimated background of less than one event. For a slant depth cut 14 kmwe 
and a muon energy threshold of 1.8 GeV, this corresponds to a muon flux of 4.01 ± 
0.50 s * a *±0.30 SJ/s xlCT 13 crn~ 2 sr -1 s _1 in the horizontal direction -0.14 < cos9 z < 
0.14. None of the 65 events have a large energy loss in the detector. We set an 
integral limit on neutrino-induced muons from AGNs and/or other sources from 
between 1.4 and 2.2 xl0~ 14 cm" 2 sr _1 s _1 depending on muon energy. 
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